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Abstract 

Nonlinear interaction of three acoustic waves in a sandy sediment is studied in the frequency 
range where there is a considerable wave velocity dispersion. The possibility of an experi¬ 
mental observation of the generation of a sound wave by two pump waves propagating at an 
angle to each other is estimated. 



1 Introduction 


Acoustic wave propagation is an effective tool for studying properties of porous media. Such 
media are known to exhibit a considerably stronger elastic nonlinearity as compared with 
homogeneous fluid or solid media [D El El HE]. This property stimulated a growing interest 
in studying nonlinear dynamics of poroelastic media (see, e.g. Refs. mmm- 

In Refs. [HI E] the generation of second-harmonic acoustic waves is studied theoreti¬ 
cally for the one-dimensional case. In the frequency range where there is no sound velocity 
dispersion, the energy and momentum conservation laws, 

001 + 002 = 00 s, ki k2 = ks (1) 

are satisfied for nonlinear interactions of three waves propagating in one direction. Such 
interactions include harmonic generation of the fundamental wave as well as the generation 
of the waves with sum and difference frequencies of the fundamental and the excited waves. 
In this case, if the dissipation is not too high, the formation of a shock wave can start, 
and a particular three wave interaction could be hidden. For more explicit observation of 
nonlinear three wave interactions it is preferable to choose the frequency range where there 
is a signihcant amount of velocity dispersion. In this case the three waves that satisfy the 
conditions ([T]) will propagate in three different directions, and no shock formation would 
occur. 


2 Theory 

In the present work the nonlinear interaction of three acoustic waves propagating in a granual 
medium at an angle to each other is studied. It has been experimentally observed that marine 
sediments exhibit a noticeable acoustic wave velocity dispersion in some frequency ranges 
[HI uni IH]. We shall use the dispersion velocity data [9] listed in Fig. 2 of Ref. [H]. The 
experimental data show a strong velocity dispersion within the frequency range between 1 
kHz and 10 kHz. We assume that at the boundary of the porous medium two waves, (cui, ki) 
and (a; 2 , k 2 ), are excited at an angle to each other. They generate the wave (cus, ks), and let 
it propagate along the x-axis. We may write 

ki cos 9i + /c 2 cos 62 = ks, (2) 

where the angles 9i and 02 are the angles between the vectors ki, k 2 and the x-axis corre¬ 
spondingly. The angles 9i and 02 depend on the choice of the frequencies ooi and 002 - We 
shall choose them in the range of the maximum velocity dispersion in order the angles 0i 
and 02 are not too small and hence don’t fall in the dissipation spreading of the waves. We 
shall be based on Fig. 2 from Ref. m and choose the frequencies, Ui = 2 tt ■ 2 ■ lO^s 
a ;2 = 27r • 3 ■ for which the sum frequency equals oos = ooi + 002 = 270 ■ b ■ 

The corresponding acoustic wave velocities are the following: ci ~ 1,58 • lO^cm/s, C 2 ~ 
1, 62 • lO^cm/s, Cs ~ 1, 72 • lO^cm/s. These experimental data yield the values of the three 
wave vectors, which allows obtaining the values of the angles 0i and 02, 0i ~ 25°, 02 ~ 18°. 

Our task is to hnd the amplitude of the wave (cus, ka) generated by the waves (cni, ki) and 
(a; 2 , k 2 ) and to estimate if the intensity of the generated wave could reach a measurable value 
at a reasonable distance starting from the fluctuation level at the boundary. To solve the 
problem we start from the continuity equations for the densities and momenta of the liquid 
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and solid phases of a sediment composed of a rigid frame and pores filled with water, see Refs. 
[ni[i3]. These equations are equivalent in the main features to the equations developed by 
Biot [Ul |T51 |TB], but they are presented in a different form with a more explicit physical 
meaning. On the basis of these equations, in Ref. [7] the equations for the densities of the 
liquid and solid phases, pf and ps, were derived (in this paper we don’t take diffraction into 
account): 
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In these equations the following variables are used. 


X = ex 


(4) 

(5) 


and the moving coordinate 


T = t — X / C. 


(In Eqs. (E]), dl]) and everywhere below the primes for x are omitted.) 
the small parameter e is introduced as 


In the relation 


( 6 ) 

o 


€ ~ v^/c ~ u^jc ~ SpilPi ~ SpJPs, (7) 

here c is the sound velocity in the sediment and f, u are the hydrodynamic velocities of the 
liquid and solid phases, 5pf, 6ps are the deviations from equilibrium values of the densities 
of the liquid and solid phases. In Eqs. ([3]), dl]) the left-hand sides are of the order of ~ e, 
the right-hand side terms are of the order of ~ e^. The introduction of the new variables 
(0, (0) actually signihes the application of the method of slowly varying wave prohle at the 
distance of the wave-length scale. 

In the equations ([3]), (HI) we wrote Pf,Ps instead of 6pf,6ps, m is the porosity; 


ks kf k} ’ 

where kf, ks and k are the bulk moduli of the fluid, mineral grains constituting the frame, 
and of the frame itself; p is the shear modulus of the frame; u = 1 — m — k/ks] P^ is the 
nonlinear part of the pressure in the fluid, = cr”^, — k/ksP^ where cx^j, is the nonlinear 
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part of the stress tensor of the frame. The nonlinear stress tensor is considered for the one¬ 
dimensional case, which corresponds to the accepted approximations. The explicit forms for 

and are given in Ref. [7] for some limiting cases. 

Let ns eliminate one of the variables, 6pf or 6ps, from the linear parts of Eqs. (jS]), (jl])) 
(let it be e.g. Sps), by snbtracting one eqnation from the other one. Note, that Eqs. ([3]), ([1]) 
allow two independent longitndinal modes, the so called fast and slow waves. As it is shown 
in Ref. [T7], the slow wave (nnlike the fast one) is a strongly attennated diffnsion mode, and 
it does not contribnte signihcantly to the sonnd held. In this single-mode approximation, 
the elimination of 6ps leads to the disappearance of the linear part of the eqnation provided c 
is the velocity of the fast wave. In the nonlinear parts the qnantity Sps is expressed throngh 
6pf with the formnla which is valid to an accnracy ~ e. 


5ps = 


OsC^G 


1 - 


m 


PfC^G^ 


^Pf- 


( 8 ) 


In this approximation we arrive at the nonlinear eqnation for an aconstic wave in a 
sediment. 
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To obtain these eqnations we took into acconnt that in sand sediments the bnlk modnlns 
ks of qnartz grains is mnch greater than that of the pore water, and in this case G can be 
evalnated as G ~ m/kf, provided m is not close to zero. 

In Eq. ([9]) the term aiDrPf that acconnts for dissipation is introdnced. Dr is the 
dissipation linear operator in the variable r which is characterized by the property 




= a{u)e''' 


( 10 ) 


where a is real and positive and it has the meaning of an amplitnde attennation coefficient if 
the coefficient Oi is taken to be eqnal to the coefficient at dpf /dx. The relation flTOj) dehnes 
the action of this operator on any fnnction of the variable r which can be represented by 
a Fonrier series or integral. An algebraic expression for Q!(a;) is a combination of physical 
parameters (complex bnlk and shear frame modnli inclnded) of a sediment, and it inclndes 
the freqnency correction fnnction introdnced by Biot |16] . 

Let ns write Eq. ([9]) in a concise form. 


ai 


dpf 

dx 


dp] 

+ OjlDrPf = 0 . 


( 11 ) 


Note the following. The coefficient 03 at the nonlinear term in this eqnation incorpo¬ 
rates contributions from the nonlinear stress-strain relations of the pore flnid, of the grains, 
constitnting the sediment frame, and of the sediment frame itself. To write Eq. CD we 
considered the nonlinear pressnre P” and stress tensor a/ir. from Eq. dH]) as expansions in 
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powers of 5p/,s, retaining qnadratic terms. The coefficient 03 inclndes all the coefficients at 
dpj/dr, and this coefficient will be estimated basing on experimental data. 

We shall use Eq. fITT]) to find the slowly varying amplitude of the generated wave ps 
(from now on we shall omit the index / at p). With the three interacting waves written 
in the form pi = \pi\ expi(kjr — Uit + ipi), i = 1,2, and ps = jpsl expi^k^x — + ps), and 

taking into account the relations ([H [ 2 ]), we can write Eq. flTTl) in the form, 


+ o|p3|e‘^’ = 0, (12) 

ax ax ai 

where a is the amplitude attenuation coefficient of the wave ps. The second term in Eq. 
ca vanishes when the phase pa attains its fixed value. The attenuation term alpsle*'^® can 
be omitted since the wave pa is generated and sustained by the waves pi and p 2 along the 
whole distance of their interaction. The attenuation of the exiting waves should be taken 
into account to estimate the real distance of the nonlinear interaction. 

One can see from Eq. 03. that the wave ps will ’’survive” and will be amplihed in case 
¥^3 = + ^2 + 7 r/ 2 , and the equation for the slowly varying amplitude ps acquires the form, 


^--‘P3|P1||P2|=0. (13) 

ax ai 

Thus, the vertex (the second-order unharmonicity), that determines the interaction under 
consideration is equal to {a 3 /ai)u 3 . From this equation we get 

Psil) - P3(0) = —UJ 3 P 1 P 2 I, (14) 

QjI 

where P 3 ( 0 ) is the amplitude (at the fluctuation level) at a; = 0 , and it may be neglected, 
since p^^l) is supposed to be much higher. 

Let us estimate the distance I at which the amplitude p 3 can reach a measurable value. 
This distance cannot exceed the dissipation lengths of the waves pi and p 2 . The amplitude 
attenuation coefficients of these waves are equal correspondingly to cki ~ 0.8 ■ 
and to 0(2 ~ 3 • (see Ref. m)- This corresponds approximately to propagation 

distances ~ 1250 cm and ~ 330 cm. This means, that I cannot exceed the distance that is a 
little more than 300 cm. 

The coefficient ( 03 / 03 ) can be estimated if the second-order nonlinear parameter usually 
denoted as B/A is known: (03/01) ~ BA~^/pc, here p is the average density of the medium. 
Water-saturated porous media are known to have a stronger nonlinearity than a homogeneous 
fluid. For example, the parameter B/A can take values ~ 8 — 12 for sandy sediments Pi], 
while in water it equals ~ 5 — 6 . In Refs. [TSl HH] the nonlinearity parameter of water- 
saturated sand was determined to be about 100. For such granual media as disordered 
packings of noncohesive elastic beads embedded in a fluid it is 10 ^—! 0 ^ times larger than in 
homogeneous fluids and solids [ 8 ]. One can also note that quadratic nonlinearity increases 
signihcantly in the presence of gas bubbles in water. Experiments IZDl showed that the 
nonlinearity parameter of water containing gas bubbles can reach 10 ^ — 10 ®. 

We shall make a numerical estimate for the distance I at which the ratio P 3 (/)/pi, 2 ( 0 ) can 
reach the value (10 ^ — 10 ^). We assume 


Pi(0) ~ P 2 ( 0 ) ~ 10 ^g/cm^, p ~ 2g/cm?, c = 1, 7 • 10®cm/s, 

ws = 27r ■ 5 • 10®s“\ B/A ~ 10. 


These data yield the distance I equal approximately to ~ 300 cm, which falls within the 
attenuation length of the exciting waves. 
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3 Summary 

Nonlinear interaction of three acoustic waves in a sandy sediment is studied . A significant 
amount of velocity dispersion in some frequency intervals allows momentum and energy 
conservation laws to be satisfied for the waves propagating not in one direction. This means 
that no shock formation would occur and hence this would not obscure a three-wave nonlinear 
process. Numerical estimates show that the wave generated by two pump waves propagating 
at an angle to each other in a sandy sediment can reach a measurable value at a distance 
realistic for an acoustic-wave experiment in a sediment. 
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